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Old models of the scattering of composite systems based on the Glauber model of multiple 
diffraction are applied to hadron-nucleus scattering. We obtain an excellent fit with only two free 
parameters to the highest energy hadron-nucleus data available. Because of the quality of the fit 
and the simplicity of the model it is argued that it should continue to be reliable up to the highest 
cosmic ray energies. Logarithmic extrapolations of p — p and p — p data are used to calculate the 
proton-air cross sections at very high energy. Finally, it is observed that if the exponential behavior 
of the p — p diffraction peak continues into the few TeV energy range it will violate partial wave 
(^-^ ■ unitarity. We propose a simple modification that will guarantee unitarity throughout the cosmic 

' ray energy region. 
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The highest energy hadronic interactions currently observable occur when a primary cosmic ray strikes an air 
nucleus. For primary protons with energy of 10 18 eV, which are easily observed in the Utah "Flys's Eye" detector, 
' the yfs exceeds that of the recently cancelled SSC. The phenomenology of the extensive air showers created by the 
hadronic cascade in the atmosphere depends critically on the hadron-nucleus and the nucleus- nucleus cross sections at 
extremely high energy. Furthermore, extracting any information about the basic hadronic interactions requires some 
model that can reliably relate the nuclear cross section to the nucleon cross section. The purpose of this paper is to 
develop such a model and check it with experimental data in the energy range in which both the hadron-nucleus and 
[ hadron-nucleon cross sections are known. 

At high energies hadron-nucleon scattering amplitudes become very simple. The amplitude is diffractive being 
i-S^ ' dominated by the imaginary part with rather weak (logarithmic) energy dependence. The momentum transfer depen- 
dence in the dominant part of the diffraction peak is a pure exponential. This is to be contrasted with the situation at 
lower energies where resonance phenomena produce rapid energy dependence and complicated angular dependence. 
In this paper we have resurrected an old model of Fishbane and Trefil, and of Franco [Q and show that this simple 
■ application of the Glauber multiple diffraction model can be used to provide reliable hadron-nucleus cross sections 
| from the basic hadron-nucleon interaction. We further argue that since no energy dependence is required in this model 
in fitting the data, i.e. the distribution of nucleons in the nucleus is independent of the momentum of the incident 
hadrons, this method should continue to be applicable to scattering at the highest cosmic ray energies. One important 
reason that the nuclear corrections are simple at high energy is that for large A the nucleus is mostly black and only 
the peripheral shell of the nucleus is sensitive to the details of the hadron-nucleon interaction; thus the nuclear cross 
sections depend only weakly on the nucleon density distributions and hadron-nucleon cross sections. It should be kept 
in mind, however, that Glauber's approximation is only valid for small scattering angles. 

In the next section we will review the Glauber multiple diffraction model for hadron-nucleus and nucleus-nucleus 
scattering, assuming a Gaussian distribution of nucleons in the nucleus. In section III, we determine the input 
parameters, the hadron-nucleon cross sections and slope parameters for p — p,p — p, n + — p, ir~ — p, K + — p and 
K~ — p. The two parameters of the nuclear model are then adjusted to fit the corresponding nuclear cross sections 
as measured by Carroll et al. 

In section IV we obtain simple fits to the energy dependence of the hadron-nucleon parameters and propose a simple 
model which satisfies partial wave unitarity and can be used to extrapolate these quantities to very high energy. These 
extrapolations are then used in section V to calculate the hadron-air and nucleus-air cross sections in the energy range 
needed in the phenomenology of extensive air showers produced by cosmic rays. In the conclusion we summarize our 
results and point out that more hadron-nucleon data are needed to improve our understanding and aid the analysis 
of cosmic ray showers. 
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II. THEORETICAL FORMULATION OF HADRON-NUCLEUS AND NUCLEUS-NUCLEUS 

SCATTERING 



Many authors including those listed in H have used the Glauber multiple diffraction approximation |2| to relate 
hadron-nucleus scatttering and nucleus-nucleus scattering to the basic hadronic interaction. In the interest of making 
this paper self contained we will review this procedure. 

The amplitude for hadron-nucleus elastic scattering can be written in the impact parameter representation as 

Fh-A(q 2 ) = ^ / d 2 be^\l - e^- A ) (2.1) 

27T J 

where Xh~A is the hadron-nuclear phase shift. The quantity Th- A — (1 — e iXh - A ) is the scattering amplitude in 
impact parameter space, often refered to as the nuclear profile function, and is bounded by unitarity because the 
imaginary part of Xh—A is positive. In this treatment we will ignore spin and treat neutrons and protons alike as 
nucleons. The essential approximation of the Glauber model is that the scattering is predominately forward and there 
is no momentum transfer to the individual nucloens. In this limit the nuclear phase shift for a particular nucleon 
configuration is given by the sum of the phase shifts for the collisions of the hadron with the individual nucleons in 
the nucleus. We believe that this should give good results for the total cross section given by the imaginary part 
of the forward amplitude and for the elastic scattering which is dominated by the forward diffraction peak. On the 
other hand, the fact that the nuclear amplitude contains no nuclear recoil makes it unlikely that this approximation 
will give good results for quasi-elastic processes in which the nucleus ends up in an excited state or fragmented. For 
cosmic ray applications there are a number of processes that fall into the general catagory of target fragmentation 
that reduce the beam momentum, but do not contribute to shower formation. Since these must be treated in the 
simulation, we see no particular reason to use the Glauber model calculation of the quasi-elastic scattering to remove 
an ill definited and perhaps incorrect fraction of these events. 

If correlations between the nucleons are ignored and we average over the position of the (A) nucleons relative to 
the center of mass of the nucleus, the overall hadron-nucleus phase shift is 

A 

e^-AW = J n&TiP^e^U »>-«( 5 -^ (2 2) 

l — l 

The contribution of each nucleon to Th-A factorizes and we obtain 

F h - A (q 2 ) = ^J <Pbe^{l -[I- J d a rT h _ N (b- ?i_)p{f)] A } (2.3) 

where p(r) is the single nucleon density, r_i_ is the component of r in the impact parameter plane and Th-N is the 
hadron-nucleon profile function which is just the Bessel transform of the hadron-nucleon scattering amplitude: 

T h - N (b) = J d 2 qe-^ b F^ N (q). (2.4) 

All of the scattering processes are at a fixed incident hadron momentum. Being able to neglect the nuclear Fermi 
momentum is another simplification that is possible at high energy. 

The multiple diffraction approximation can also be applied to nucleus- nucleus scattering || (see also Refs. Q and 
). In this case the overall phase shift is, according to the Glauber assumption, simply the sum of the phase shifts 
of A x A' individual nucleon- nucleon interactions. In this case we must average over the positions of the nucleons 
in both nuclei. Unlike the hadron-nucleus case, however, the integrals do not factorize and lead to very complicated 
expressions even for the 2x3 case. It can be shown || that in the optical limit, that is both A, A' 3> 1, many terms 
can be neglected and the elastic scattering amplitude can be approximated by: 

F A ,_ A (q 2 ) = ^-f d 2 be^ b {l - [1 - J d 3 r J dV p A (r)T N - N (b + f x - 7 ± )p A , (f (2.5) 

In practice the optical limit is a good approximation for A, A' > 10 M. From this expression, using the optical 
theorem, we otain the total cross section 

<ta-a> = T ImF A'-A{0) (2.6) 
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and the elastic cross section 



a A-A 



= J dQ k \ F AI _ A {q 2 ) | 2 . (2.7) 



In this work we will use a Gaussian single nucleon density 

Mr1 = ^?e"W. (2.8) 

This is of course only valid for A > 1 and for A = 1 the density should be a <5-function. While other distributions 
such as Saxon- Wood are probably more realistic for large A and shell wavefunctions are more appropriate for small 
A, this form seems to work very well. When the density given in Eq. 2.8 is combined with the observed form for the 
hadron-nucleon scattering amplitudes, Eqs. 2.6 and 2.7 can be evaluated analytically. We assume that R(A) has the 
following form: 

R(A) = RqA' 1 . (2.9) 

If we use experiment to determine T^-n, the only free parameters in this model are Rq and 7. These could of course 
depend both on energy as well as particle type, although the simplest picture is one in which the density only reflects 
the average positions of the nucleons in the nucleus. 

The hadron-nucleon elastic amplitude at small momentum transfer are well fit by a Gaussian form. More compli- 
cated dependence appears at larger angles, but it is the simple small angle region that makes the dominant contribution 
to the elastic cross sections. This remains the case up to the highest experimental energies available at the Tevatron 
collider. The form of the scattering amplitude is then 

F h . N = k -^{p + i)e^ (2.10) 

where B is the slope parameter (t = —q 2 ) and p is the ratio of the real to imaginary part of the amplitude. Here we 
have assumed that for small q 2 , p is a constant. At high energy p is quite small, typically 0.10-0.14, and we take it to 
be zero in our work. The resulting hadron-nuclcus absorptive cross section is 

and the total cross section is 

,r_ A = 2 ,f° dx (l±^± = 20 £(A^ ( , 12) 



n=l 



where 



and 



/3 = -n{2B + R 2 {A)) (2.13) 



a 



T 

h-N 



I>=-|f. (2-14) 

Finally, we include the prediction for the quasi-elastic scattering. While we are doubtful about its reliability, it turns 
out be be small and we will be force to use it in our fit to accelerator data. 
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.„_ , - ■ . , dx- £±|> ofc A (2.15) 



where 7 = ^ 



The nucleus-nucleus cross sections can be obtained from Eqs. 2.11 and 2.12 by letting A — ► AA' and using 

13 = -tt(2B + R{A) 2 + R(A') 2 ) (2.16) 
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and 



D=-yL (2 . 17) 



The following identity Q 



r D (l + x) M -l 

Lim M ^oo / dx K ' = Ei(MD) - ln(-MD) - C (2.18) 

Jo x 

where C is Euler's constant, can be used to evaluate these cross sections for large A or AA'. Up to this point we 
have neglected position correlations of the nucleons in the nucleus. Indeed, the nucleon coordinates f are subject 
to the center of mass constraint -4 53i=i ^ = ^cm- It can be easily shown, by transforming to the cm coordinates 
f? m = f — -j YLf=i that the effect of the cm constraint on the elastic scattering amplitude, for the case of Gaussian 
nucleon density distributions, is simply 

FZ m A (q)=ei 2! ^ 1 F^ A ( q ) (2.19) 



*A-A' 



(q) 



2, R{A) 2 R(A') 2 , 

4A> > Fa- 



(2.20) 



Thus the total cross sections are unaffected by the cm constraint. The inelastic and elastic cross sections however are 
modified, although as it is shown in Fig.l the correction is at most a few mb and decreases as ~ A~ x l z . We find 



Acr p-A 



2 4±± 

jV=l M=l 



and for A > 1 
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2B+R(A) 2 
R.{A) 2 



2MN 1 



(2.21) 
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(2.22) 



A recent model developed by Engel et al. (hereafter the EGLS model) [Q and discussed in more detail by Fletcher 
et al. |§ is based in part on the work of Durand and Pi |J, which is a semi-classical treatment motivated by the 
parton model. The EGLS model is a probabilistic approach to hadron-nucleus and nucleus- nucleus inelastic scattering, 
treating the individual hadron-nucleon interactions incoherently, and summing probabilities of interaction instead of 
amplitudes. While this treatment is somewhat different from ours, for the special case of purely imaginary nuclear 
phase shift the results are identical. 



III. CROSS SECTION FITS 

The input parameters of Glauber's model are the hadron-nucleon cross sections uf^_ N and nuclear slopes Bh-N- 
We determine these at energies of 60, 200, and 280 GeV. The total cross sections are found in a compilation of data 
||. Some sample fits are, (s is in GeV 2 and a T is in mb) 

(Tj + _ p = 23.1 + .525(ln S /112) 2 (3.1) 

and 

CTj-_ p = 24.0+.582(lns/170) 2 (3.2) 

with a \ 2 = ifj- and \ 2 — ^rf correspondingly. There are some values of the slope parameters in the literature 
fiol |l2| but we obtain most values by fitting directly the elastic differential cross section data (at low momentum 
transfer) with an exponential form: 

d ^gL=Ae Bt . (3.3) 
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Sample fits of these values are, (s is in GeV 2 and B is in GeV 2 ) 

B„-_ p = 9.53 - 0.7051ns + 0.136(lns) 2 (3.4) 



B w +_ p = 7.90- 0.1551ns + 0.076(lns) 2 . (3.5) 

Ayres et al. jl3| have measured ir — p and K — p differential elastic cross sections at 50, 70, 100, 140, 175, and 200 GeV. 
There are some data for tt~ — p at 360 GeV [H and for n — p at 200 GeV jl5|. By interpolating and extrapolating 
in the case of K — p we obtain Bh- P at the required energies. The only free parameters in our model are Rq and 
7 that appear in the nucleon density distribution. The most recent hadron- nucleus data is due to Carroll et al. [0. 
The comparison of our model with their data is complicated by the fact that they removed the quasi-elastic cross 
section from their absorption cross section. While this is a small correction to the absorption data it was done in the 
process of extrapolating the data to zero momentum transfer and hence we cannot obtain the total absorption cross 
sections from their published data. As a result, we have two courses of comparison. First, we can simply ignore these 
small corrections and fit their data with our model. The alternative is to use the Glauber model to calculated the 
quasi-elastic cross sections. This has the disadvantage of subtracting a small but unreliable term from the quantities 
that we believe can be calculated accurately. 

We first ignore the quasi-elastic cross sections and simply fit the 108 data points of Carroll et al. H for scattering 
of p,p, K + , K~, 7r + , 7r~ off Li, C, Al, Cu, Sn, and Pb targets at 60, 200, and 280 GeV incident hadron energy with 
our two free parameters. These results are shown in Figs. 2 and 3. We obtain, 

R = 3.89 GeV- 1 (3.6) 



7 = 0.31 (3.7) 

with x 2 /d.f- — 0.25. This x 2 was calculated using the quoted error of the data which includes an estimated systematic 
error of 3%. The small value of our \ 2 is an indication that our fit is not sensitive to the type of systematic error that 
exists in the data. If we use only the statistical error of 1% as quoted by Carroll et al. j| to calculate \ 2 we obtain a 
X 2 ~ 2 per degree of freedom. 

If we use the Glauber model to remove the quasi-elastic scattering from our fit, we obtain, 

R = 4.74 GeV^ 1 (3.8) 



7 = 0.28 (3.9) 

with x 2 /d.f. = 0.34. While this is not as good a fit as the one described above, it is quite acceptable. Note that there 
is a noticeable change in the parameters. 



IV. UNITARITY CONSTRAINT AT HIGH ENERGY 



The data for the total nucleon- nucleon cross section a T and nuclear slope B up to Tevatron energies 15 can be fit 
with quadratic logarithmic forms and at infinite energy they grow as ln(s) 2 . Our fits, shown in Figs. 4 and 5 are: (s 
is in GeV 2 , a T in mb and B in GeV -2 ) 

cr T „ p = 38.46 + 0.41(ln(s/118)) 2 (4.1) 

and 

B p - p = 9.25 + 0.371n(s) + 0.01(ln(s)) 2 (4.2) 

T 

where we have assumed that at p — p and p — p are identical at collider energies. Note that for s — > oo, T(0) = — > 
8.35. In fact, the unitarity limit is saturated at y/s ~ 2.5 TeV (see Fig. 6). To prevent unitarity violation we suggest 
that at high energies the elastic scattering amplitude deviates from the pure exponential form in t that is commonly 
used at small scattering angles. We set, 

f(t,a) = £o T (8)e S # i J [a(8)V=i\ (4.3) 
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Since dJo(x)/dx\ x =o — 0, B is still denned as the nuclear slope of the elastic amplitude. The profile function that 
follows from Eq. 4.3 becomes smeared with a modified Bessel function: (see Fig. 6 ) 

a 2 +fa 2 Oh 

and a(s) is chosen such that T(b, s) < 1 or a > [2B m(^rg )] 1 ^ 2 - For ^/s < 2.5 TeV we set a — so that our cross 
section fits are not modified since they are well below this energy. At higher energies the above profile function leads 
to 

a T h _ A = -2/3 / dx{l - [1 + D'e-*I (\V^)} A } (4.5) 
Jo 



POO 

oT-A = ~P / dx{l - [1 + D'e-*I Q (\^)} 2A } (4.6) 
Jo 

where (3 = -tt{2B + R(A) 2 ), A = 2a y/n/ - /3 and D' = De x / 4 . The nucleus-nucleus case can be obtained as before, 
in the optical limit, by substituting A -> AA' and (3 = -ir(2B + R(A) 2 + R(A') 2 ). As it is shown in Fig. 6 for 
A = Air, the p — A profile function remains unchanged for small impact parameter, (saturating the unitarity limit 
for large ^4) but increases at large b, due to the smearing of the nucleon-nucleon profile function. This results in an 
increase of the total and inelastic p — A cross sections, as shown in Fig. 7. In these plots we have used the minimum 
value of a(s) that satisfies unitarity. 

Our model of course is by no means unique since there may be other ways to preserve unitarity. The choice of the 
Bessel function is not essential and other modifications such as exponentials etc. should produce the same qualitative 
results. 



V. EXTRAPOLATION TO COSMIC RAY ENERGIES 



The development of cosmic ray showers in the atmosphere depends strongly on the size of the relevant hadron- 
nucleus and nucleus-nucleus cross sections at very high energies E ~ 10 18 eV, y/s ~ 40 TeV. Both the interaction 
length in the atmosphere and the number of nuclcons in the nucleus that participate in a collision are a function of 
the absorption cross sections. 

In the previous sections we have shown that the hadron-nucleus cross sections can be calculated reliably from the 
basic hadron-nucleon interaction. Since Glauber's analysis does not depend explicitly on the energy scale at which 
the scattering takes place we expect that it is still applicable at very high energy. On the other hand, only for p — p 
does experimental data approach modest cosmic ray energies. Even this data at y/s — 1.8 TeV (Ei a b — 1.7 x 10 6 GeV) 
must be extrapolated over three decades in energy to reach the highest cosmic ray energies of E p > 10 9 GeV. Because 
of the observed logarithmic behavior of the cross section and slope parameter, this procedure should be fairly reliable 
and data that will eventually come from the LHC should provide direct input for all but the very highest energies. 
The situation for the K — N and tt^ 1 — N interactions is much less satisfactory. Data for these reactions exist only 
up to a few hundred GeV lab energy with little hope of future measurements above the TeV range even with the 
LHC. While external beam experiments up to the Tevatron's maximum energy would improve our understanding of 
these reactions, it is unlikely that energy extrapolations over 7 or 8 orders of magnitude can be trusted, which means 
that the forms given in Eqs. 3.1 through 3.4 cannot be used at cosmic ray energies. It appears that the best method 
for determining — N and — N scattering at cosmic-ray energies is to use the independent quark model result 

o 

a 7r-N = a K-N = ^ a N-N (5-1) 

Although this is only approximately true at low energy, it is expected to improve at very high energy. The hadronic 
nuclear slopes are taken equal. These approximations are probably adequate considering the fact that the hadron- 
nucleus cross-sections are relatively insensitive to the precise value of -B/j-jv and as it is shown in Figs. 8 and 
9 for p — Air. In Fig. 7 we plot our prediction for the p — Air cross section as a function of center of mass energy 
and compare with a fit by Block and Cahn of [|l7j. Based on theoretical arguments they assumed that the total 
proton-proton cross section has the form 

a T p _ p = A + (3 (ln 2 (j-) - + csin(^)^- 1 + ^cos^)^ 1 (5.2) 



G 



The parameters from their fit are A = 41.74 mb, = 0.66 mb, s = 338.5 GeV 2 , c = 0, D = -39.37 mbGeV^ 1 ""), 
and a = 0.48. The nucleus-nucleus cross sections can also be calculated at cosmic ray energies based on the nucleon- 
nucleon interaction using Eqs. 2.11 and 2.12 with the appropriate substitutions. However, there are not enough 
experimental data on nucleus-nucleus cross sections to check whether the Glauber model works as well as in the 
hadron-nucleus case, especially since for practical purposes one is forced to use the optical limit approximation. 



VI. CONCLUSION 



We have shown that Glauber's model of multiple diffraction accurately predicts hadron-nucleus cross sections from 
the basic underlying hadron-nucleon interaction. We have obtained excellent agreement with the experimental data 
for a variety of incident hadrons and target nuclei in the energy range of 60-280 GeV with just two free parameters. 
External beam experiments at 800 GeV on nuclear targets would provide a useful check of the validity of these 
methods over a much broader range of energies. Since Glauber's analysis does not depend explicitly on the energy 
scale of the scattering process and it is unlikely that the distribution of nucleons inside a nucleus would change with 
energy, we expected this method to continue to be applicable at very high energies. The p — p data extend to a large 
energy range and can be extrapolated to cosmic ray energies. Using Glauber's model we can then obtain the relevant 
nucleon-nucleus cross sections at these energies. We can use the independent quark model to estimate the ir — N and 
K — N cross sections, but have not yet developed a model for the nuclear slope. 

Finally, we would like to emphasize that extrapolating our fits to available data on and B p - p to energy y/s 
greater than a few TeV violates unitarity. This means that above the current Tevatron energies and the energy 
available at the LHC we should expect some change in the shape of the nucleon diffraction peak and its related profile 
function. This is to be contrasted with the behavior at lower energies where the Gaussian profile function height and 
radius have grown logarithmically with s. 

We propose a simple model in which the p-p elastic scattering amplitude deviates from a pure exponential in the 
momentum transfer and satisfies unitarity. The resulting absorption cross sections are increased appreciably at cosmic 
ray energies. How nature actually chooses to satisfy unitarity awaits experiments in the few TeV energy range. 
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FIG. 1. Correction to the p— A inelastic cross-section due to the center of mass constraint as a function of A at ^/s — 10 2 GeV. 
The dotted line is the large A approximation. 

FIG. 2. Fits to p — A, tt + — A, and K + — A absorption cross sections at incident hadron momenta of 60, 200, and 280 GeV/c. 
Shown are the data by Carroll et al. and our fitted values connected by straight line segments. 

FIG. 3. Fits to p — A, n~ — A, and K~ — A absorption cross sections at incident hadron momenta of 60, 200, and 280 GeV/c. 
Shown are the data by Carroll et al. and our fitted values connected by straight line segments. 

FIG. 4. The p — p total cross section as a function of s. The solid line is our fit of Eq.4.1, and included are the high energy 
p — p collider data which we assume are equal to p — p at these energies. 

FIG. 5. The p — p nuclear slope B as a function of s. The solid line is our fit of Eq.4.2, and included are the high energy 
p — p collider data which we assume are equal to p — p at these energies. 

FIG. 6. Profile functions for p — - p and p — Air scattering at s = 10 9 GeV 2 . Shown are r p _ p the unitarity violating profile, 
the Bessel-modified profile Fp_ p and the p — Air profiles that follow from them using Glauber's approximation. 

FIG. 7. Inelastic p — Air cross sections as a function of cm energy. The solid line is the Bessel-modified cross section which 
is higher than the unitarity violating cross section above ~ 2.5 TeV (dashed segment). Also shown for comparison as a dashed 
line, is the Block and Cahn parametrization of Eq. 5.2 

FIG. 8. Dependence of the p — Air inelastic cross section on the nuclear slope of the proton-proton elastic scattering, keeping 
the p — p cross section constant at 40 mb. 

FIG. 9. Dependence of the p — Air inelastic cross section on the proton-proton cross section, keeping the p — p nuclear slope 
constant at 10 and 20 GeV~ 2 . 
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